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The heat capacity of some ferrimagnets has additional structures like a 
shoulder in the Schottky-like peak, or emergence of a second peak when an 
external magnetic field is applied. It is shown here that as long as spin 
wave-spin wave interactions are ignored in a ferrimagnet, the ferromagnetic 
and antiferromagnetic elementary excitation spectra give rise to two inde- 
pendent heat capacity peaks, one enveloped by the other, which add up to 
give the peak for the total system. Taking this into account helps under- 
stand the additional structures in the peaks. Moreover, the classification of 
ferrimagnets into predominantly antiferromagnetic, ferromagnetic, or a mix- 
ture of the two is shown to be validated by studying them under additional 
influences like dimerization and frustration. Because these two are shown 
to influence the ferromagnetic and antiferromagnetic dispersion relations - 
and hence the quantities like heat capacity and magnetic susceptibility - by 
different amounts, the characterisation of ferrimagnetic systems (1, |), (|, 1) 
and (§, §) is brought out more clearly. Both these influences enhance anti- 
ferromagnetic character. 

PACS numbers: 75.10.Jm, 75.50.Ee 
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1 Introduction: 



Thermal properties of ferrimagnetic chains were theoretically investigated re- 
cently employing various methods like the modified spin- wave theory (MSWT) |], 
0, density matrix renormalization group (DMRG)0, quantum Monte Carlo 
method (QMC)|l], [|, and Schwinger boson mead field theory (SB)[[|. The 
specific heat C v and magnetic susceptibility \ were shown to depend upon 
temperature as T 1 / 2 and T~ 2 respectively at low temperatures, and yT and 
C v were shown to have, respectively, a rounded minimum ||, |[] and a 
Schottky-like peak|| |J at intermediate temperatures. The spin correlation 
length was shown to have a Independence at low temperature Q, |5|. The 
modified spin wave theory, modified either by including Takahashi constraint 
or by including higher order corrections in the spin wave theory, was also 
shown to give results in surprisingly good agreement with those from quan- 
tum Monte Carlo method in the thermodynamic limit for this system]!], [|. 

Ferrimagnetic systems have been classified into three categories: one with 
a predominantly Ferromagnetic (F) character, the second with a predomi- 
nantly antiferromagnetic (AF) character, and the third with a mixture of the 
two. If the two spins constituting a ferrimagnet are Si and s 2 with s\ > s 2 , 
then it is conjectured || that the systems in the first category are those 
with si > 2s2, those in the second category have si < 2s2, and the systems in 
the third category have s\ = 2s2- We call it Yamamoto classification. Thus 
the system (§, |) is ferromagnetic in character, is antiferromagnetic, 

and (1, |) has a mixture of the two characters. 

The xT vs temperature curves of different ferrimagnetic systems look 
alike - a rapidly decaying ferromagnetic part at low T, a rounded minimum 
at intermediate T and a linearly increasing antiferromagnetic part at high T 
- except that some systems like the predominately ferromagnetic (|, |) have 
a smaller rate of increase with temperature after the minimum compared to 
others. The heat capacities of the three systems also have qualitatively the 
same shape: a T 1 / 2 dependence at low T, a Schottky-like peak and a decay 
at large temperatures. The MSW results on (|, |) system, however, show a 
shoulder in the heat capacity below the peak temperature, which has been 
explained as being a result of the deficiency of the theory 0]. When an exter- 
nal magnetic field is applied, a second peak appears at low temperatures ||]. 
It is, therefore, quite likely that the heat capacity of ferrimagnets has, under 
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suitable conditions, inherent structures like shoulders and double peaks. We 
would like in this paper to create such conditions and understand the nature 
of the heat capacity peak of ferrimagnets. 

In a recent workQ, we used a zero-temperature linear spin wave theory 
to study ferrimagnetic systems under the effects of dimerization and frus- 
tration, represented by the parameters S and a. Dimerization as well as the 
frustration were shown to affect the dispersion curves by either pushing them 
up or pulling them down in energy. Since this would change the gap that 
determines the heat capacity peak, these two effects are expected to influence 
the peak structure. It is also possible that they may give rise to additional 
peaks in heat capacity. 

There were two distinct values of the frustration parameter; the critical 
value a c = 2 ( S i+s 2 ) a ^ wmcn the l° n S ran g e order is destroyed, and point 
oc* = ^ttttt, that heralds transition from a commensurate ferrimagnetic 
state to a spiral state. The heat capacity and susceptibility are also expected 
to show the tell-tale signs of the transition induced by a. 

Linear spin wave theory is known to give a fair picture in the case of 
ferromagnetic chains and gives only a qualitative picture for antiferromag- 
netic chains @. It has also been shown that it gives sufficiently good results 



for a ferrimagnetic system ||T0|, |TT| , [12], [13| . In a frustrated system, it has al- 
ready been argued that the LSW theory yields satisfactory results at least in 
the limit of small frustration [O. The use of LSW for larger values of frus- 
tration is indeed unreliable, but it is expected to give a qualitative picture 
that we are seeking here. As noted above, there was no remarkable improve- 
ment in the results at non-zero temperatures when linear spin wave theory 
was modified either by introducing a Takahashi constraint]^] or by including 
higher order corrections in the spin wave theory]!], 0. Nor were the results 
at zero temperature any more improved by considering non-linear spin wave 
theory |T^]. The LSW theory is therefore expected to be valid in obtaining 
qualitative results even in the presence of dimerization and frustration in the 
thermodynamic limit. 

In this paper we will study alternating spin systems formed with spin 
values (si, S2) using linear spin wave theory. We investigate the temperature 
dependence of magnetization, specific heat, free energy, and susceptibility. 
We would like to see how the predominantly ferromagnetic, antiferromagnetic 
or mixed characters of the systems are brought out by subjecting them to 
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dimerization and frustration. In particular, we would like to see the effect 
of the frustration-induced ferrimagnetic-to-spiral state phase transition on 
heat capacity and magnetic susceptibility. It is suggested that such peculiar 
effects will help identify the presence or absence of spin-Peierls dimerization 
and frustration in real low-dimensional ferrimagnetic systems. 



2 Linear spin wave theory: 

We consider a chain with spins s% and s 2 ( s\ > s 2 ) sitting on alternating 
sites with the possibility of lattice distortion leading to dimerization, and of 
competing antiferromagnetic nearest and next nearest neighbor couplings, J\ 
and J 2 respectively. A two sublattice model of this system may be described 
by the Hamiltonian 

h — X Jl^f ■ ^h-i + ^2 X [sf ' s(+2 + sf +l ■ Sf +3 , (1) 

i i 

with S A belonging to one sublattice containing spins s\ and S B to the second 
sublattice containing spins s 2 . The interaction J\ = Ji(l+ ( — describes 
the alternate weaker and stronger couplings between two adjacent sites that 
may come about because of a spin-Peierls dimerization of the lattice. 5 is 
the diemrization parameter. The degree of frustration is given by the ratio 
a = 

■Ji 

In the standard linear, non-interacting spin wave analysis, the above 
Hamiltonian is written in terms of bosonic spin-deviation operators with the 
help of Holstein-Primakoff transformations, linearized and then diagonalized 
in terms of normal mode operators to 

H = e g + J2 [Ei{k)a\a k + E 2 (k)Pt(3 k ] (2) 

k 

with Eg is the ground state energy, E lt2 are the acoustic and optical mode 
energies and the sum runs over half the Brillouin zone. The nomenclature 
used here is the same as in our previous report 0. 
The free energy of the system is 

F = e g + J2[n a Ei(k)+n b E 2 (k)}, (3) 
k 
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where n a and n& are the bose distribution functions for the two modes. It was 
reported earlier that the free energy decreases with increasing temperature || 
in a ferrimagnetic chain. 

The static susceptibility x is 




where M = J2{^a + &b)- The product \T shows a minimum at intermediate 

i 

temperatures, indicating that both the ferromagnetic and antiferromagnetic 
modes coexist in ferrimagnets. 

The sublattice magnetization is calculated by taking the thermal average 
of the occupation number as 

Mi = sx-£n£ (6) 

k 

M 2 = J2n b k -s 2 (7) 

k 

and 

M tot = M 1 + M 2 = s 1 -s 2 -Y,( n a- n h ) (8) 

k 

where fi% = vf,{n a + n b + 1) + n a , n b k = vl(n a + n b + 1) + n b and v\ is the 
coefficient in the Bogoliubov transformation; see Ref.|§. 

Specific heat C v has a Schottky-like peak. In a text-book two level system, 
the peak is a result of the gap, and an increase in this gap makes the peak 
shift to higher values of temperature, and a decrease makes it shift to lower 
values. 

A ferrimagnet has two low lying excitations, one ferromagnetic in char- 
acter (the acoustic dispersion) and the other antiferromagnetic (the optic 
dispersion curve). The Schottky like peak in the heat capacity is understood 
to be due to the two dispersion curves. In principle, because of the two ex- 
citation dispersion curves, one would expect this to behave like a three level 
system. But the heat capacity of a ferrimagnet does not follow the pattern of 
a 3 level system. Instead, from the non-interacting spin wave theory results 
of ferrimagnets, the contributions of the ferromagnetic and antiferromagnetic 



5 



dispersions to the heat capacity of a ferrimagnet seem to be additive. The 
heat capacities calculated for the modes separately has each a Schottky-like 
peak structure, and the two add up to give the heat capacity of the ferrimag- 
net. This is also expected from the fact that in the absence of spin wave-spin 
wave interactions, the free energy is a sum of the ferromagnetic (F) and anti- 
ferromagnetic (AF) normal mode energies, as in Eq.(Q). The F and AF heat 
capacity peaks, which have the usual T 1 / 2 and T 2 dependence, respectively, 
at low temperatures, are mostly overlapping due to which the resultant heat 
capacity has a single peak. The additive nature of the two contributions 
is also reflected in the dependence of susceptibility on temperature. With 
increasing temperature, the ferromagnetic contribution decays earlier, while 
the antiferromagnetic one increases and persists up to higher temperatures 

1 1, i I 

The behavior of these physical quantities with temperature in the pres- 
ence of dimerization 6 and frustration a parameters for different system 
(sij s 2); (1, 5)) (§j §) an d (§, l)i shall be discussed below. The effect of dimer- 
ization alone is discussed in section III, and that of frustrating alone in section 
IV. 



3 Dimerized chains 

The free energy is known to decrease with temperature. We find that it 
decreases with dimerization also for all (s\,s-2) systems, as shown for the 
(1, |) system in Fig.(l). The LSW results show that the free energy scales 
as T 1 / 4 , and 5 as 5 2 . 

Total magnetization decreases with increasing T and 5, with a slower 
decrease for higher dimerization. This decrease is illustrated in Fig. (2), and 
scales as T , and 5 as n^r- 

Fig. (3) shows the effect of dimerization on the heat capacity of the (1,|) 
ferrimagnetic chain which is supposed to have a mixed F and AF characters. 
Fig. (3a) gives the heat capacity without dimerization, and Fig. (3b) with 
5 = 0.9. The figures show the heat capacities when the contributions of the 
ferromagnetic and antiferromagnetic excitation branches are taken separately 
and also when taken together. As expected, the total heat capacity is the sum 
of the heat capacities from the individual excitation branches in each case. 
It comes out that while both the F and AF peaks shift to lower temperatures 
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with dimerization, the effect is different on the two. The F peak shifts by a 
larger amount than the AF peak, and the magnitude of the AF peak increases 
with 5, while that of the F peak decreases. This is true for the other two 
ferrimagnetic systems (1,1) and (|, ~) also, and appears to be a universal 
feature of a ferrimagnet, be it predominantly F, AF or a mixture the two 
characters. The increase and decrease in the AF and F peak magnitudes 
clearly shows that dimerization increases the antiferromagnetic character and 
decreases the ferromagnetic one. The net effect of this different effect is the 
emergence of the second peak in the total heat capacity, and a decrease in 
the magnitude of the main peak for the (1,|) system, as shown in the inset. 
The predominantly F system (§, |) has this effect so accentuated that the 
shoulder, which in the LSW and MSW analyses is present even at 5 = 0, 
turns into a double peak, as shown in Fig. (4a). The predominantly AF 
system (1,1), on the other hand, shows no second peak or even a shoulder, 
as in Fig. (4b). It is interesting to note that the difference arises despite the 
fact that the F and AF dispersion curves are lowered by an equal amount by 
dimerization. It appears to be a result of a more involved interplay of the 
^-dependence of the excitation energies and the bosonic numbers in the free 
energy. 

The shift of the peaks to lower temperatures can be understood in terms 
of the energies of excitation. The heat capacity peaks are indeed Schottky- 
like peaks in that they are a result of the gaps in excitation spectra. The 
Schottky peak in the heat capacity of a simple text-book two-level system 
with a gap A between the two levels has a position that shifts with the 
magnitude of A. As A decreases, it shifts to lower temperatures, and vice 
versa. For a ferrimagnetic chain, contribution to the peak comes from all the 
modes which have well-defined dispersions in k-space, each mode having a 
different gap. With dimerization, the dispersion curves are lowered in energy, 
giving rise to a shift of the peak to lower temperatures. 

That the shifting is a result of the changes in the excitation spectra is also 
supported when the effect of dimerization on heat capacity is calculated by 
using a different parametrization of dimerization. We had earlier introduced 



as another possibility^, [T^] . The reasons are given in that reference. This 



parametrization has an opposite effect on the dispersion curves. Instead of 
decreasing in energy with 6 as in the case of the usual J(l±<5), the dispersion 
curves are now pushed up in energy with increasing 5, with a nonlinear 
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dependence on 5. Fig. (5) shows the heat capacity of the (1,|) chain with 
this parametrization. The Schottky peaks now shift to higher temperatures 
with increasing 5. Also, the individual F and AF peaks respond differently 
to 5. Just as the rate of increase of the AF dispersion curve with S is much 
larger at higher 5 values than that of the F dispersions, the AF Schottky 
peak shifts by a much larger amount than the F peak. Consequently, there 
comes a time when the AF peak is pushed so far to the high temperature 
side that a double peak structure develops. This appears more clearly in 
the case of the (|, |) ferrimagnetic chain in which ferromagnetic character 
dominates. This is shown in Fig. (6). The AF peak is so strongly pushed 
to higher temperatures that the F and AF peaks almost separate from each 
other. Correspondingly, in the dispersion curves of this system, the optic 
mode is pushed up by a far larger magnitude than the acoustic one. This 
establishes then that the heat capacity peaks reflect the effects an external 
or internal influence has on the dispersion curves. 

The change of character upon dimerization shows up in the susceptibility 
of ferrimagnets also. Fig. (7) shows %T vs T, the low temperature part of 
which is dominated by the F contribution and the high temperature part by 
the AF contribution. The minimum at intermediate temperatures is a result 
of a comparable contributions of the two. With dimerization using J(l ± 5), 
the increase in the slope at high temperatures testifies to the increasing AF 
character of a ferrimagnet. This increase is the largest for the predominantly 
AF (1,1) system; Fig. (7c), smaller for the (1,|) system; Fig.(7a) and the 
smallest for the predominantly ferromagnetic (|, |) system; Fig. (7b), showing 
again the validity of the Yamamoto classification. 



4 Frustrated ferrimagnetic chains 

The effect of frustration in ferrimagnetic systems has been investigated [E 



T2|| at zero temperature. Two key values of the frustration parameter were 
identified; namely, a c which marks complete destruction of the long range 
order in the system and a* which marks the onset of phase transition to 
a spiral spin structure. At any finite temperature, we find that within the 
linear spin wave theory the two are temperature independent. We will study 
the thermal effect of a frustrated chain in the absence of dimerization. 

Free energy decreases with frustration. This decrease is faster for higher 
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a. It was reported earlier [|, [L2| that at T = the ground states energy- 
increases with frustration and then starts to decrease for values of frustration 
close to a c . This behavior is illustrated in Fig. (8) for the spin system (1, ~). 

The effect of frustration on heat capacity will again be discussed for the 
F and AF contributions separately. Frustration causes the AF C v peak to 
slightly increase in size and shift to lower temperatures. But its effect on the 
C v peak due to the F excitation modes is more dramatic. Here the transition 
to the spiral phase at a* has its impact. Before a* the peak decreases in 
size and shifts to lower temperatures. Beyond this value, when the system 
goes into the spiral phase, a further increase in frustration causes the peak 
to increase in size and shift to higher temperatures. All of this is a result of 
the mode-softening changes that the F dispersion curves go through, as seen 
in Fig. (9). The total C v being a sum of the AF and F contributions in the 
non-interacting spin- wave theory changes with a as shown in Fig. (10). 

In the case of the (1, |) system, the a = chain having a mixed F and AF 
characters, the C v curve has an AF-like T 2 dependence at low temperatures. 
With added frustration, the F peak separates out at low temperatures making 
the total C v curve have an F-like T 1 / 2 dependence at low T. The F peak 
diminishes in size with increasing a, almost vanishing at a*, beyond which it 
regrows and eventually merges with the AF peak. At this point the C v curve 
again has an AF-like T 2 dependence at low temperatures. This is shown in 
Fig.(lOa) 

In the case of the predominantly F system (§, |), the heat capacity of 
which has a low temperature shoulder due to the F peak even in the absence 
of frustration, the changes affecting the F peak are very prominent, as in 
Fig. (10b). As a increases, the F peak separates out and diminish in size 
before a*, after which it rises and eventually merges with the AF peak. The 
net C v also experiences a low-temperature cross-over from a ferromagnetic 
T 1 / 2 dependence at zero frustration to the antiferromagnetic T 2 dependence 
at a beyond a* in the spiral state. 

In the case of the predominantly AF system (§, 1), the F peak in itself ex- 
periences the same changes as in the (1, §) system, but it remains submerged 
in the AF peak, and the net C v never shows a second peak, as in Fig. (10c). 
It, however, also experiences the cross-over from T 2 to T 1 / 2 behavior at low 
temperatures. 

The effect of frustration on x is the same as that of dimerization. The 
minimum shifts to lower temperatures and the slope of the high temperature 
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tail increases as shown in Fig. (11). 

In short, the effect of dimerization and frustration on the thermal behavior 
provides a nice evidence of the validity of the Yamamoto classification of 
ferrimagnetic chains. 

We thus conclude that the heat hapacity of a ferrimagnet can have a 
more complicated structure than a simple Schottky-like peak. It can have a 
low temperature shoulder or a double peak structure. The peak can shift to 
lower or higher temperatures and can also change in size depending upon the 
influences the system is undergoing. All of this can be understood in terms of 
the component ferromagnetic and antiferromagnetic peaks that retain their 
individuality as long as the spin wave-soin wave interactions are ignored. 
The Yamamoto classification of the ferrimagnetic systems into predominantly 
ferromagnetic (si > 2s 2 ), antiferromagnetic (s 1 < 2s 2 ) and ferrimagnetic 
(si = 2s 2 ) appears to have a basis in this characteristic. 

A. A. would like to thank Wu Congjun for discussion and help. 
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Figure captions 

Figure 1: Dependence of the free energy on T for several values of the 
dimerization parameter 5 in one-dimensional spin system (1,|). A similar 
behavior was found for the other two systems (|, |) and (|, 1). 

Figure 2: The total magnetization of the alternating spin chain (1,|) vs 
temperature for different values of the dimerization parameter 5. A system- 
atic behavior was observed for the other two systems (|, \) and (|, 1). 

Figure 3: Heat capacity for the chain (1,|) as a function of T . The heat 
capacities calculated with only the ferromagnetic (acoustic) branch of the 
elementary excitation spectrum and with only the antiferromagnetic (optic) 
branch are shown separately. Their sum is also shown which is exactly the 
heat capacity when the two excitation branches are taken into account to- 
gether. Figure (a) is for 5 = 0.0, and (b) for 5 = 0.9. The inset for total 
C v for the two values of 5 highlights how the Schottky-like peak shifts in 
temperature, reduces in size, and acquires a shoulder. The curves are shown 
for the system (1,|). 

Figure 4: (a) Heat capacity of the ferrimagnetic chain (|, |). This being 
a predominantly F system, the Schottky-like peak has a shoulder even in 
the absence of dimerization. The effect of dimerization is strong enough to 
separate out the F and AF peaks to result in a two-peak structure. Fig.(b) 
is for the predominantly AF system (|, 1) which fails to show a shoulder or 
a double peak structure even at large dimerization. 

Figure 5: The dependence of heat capacity on T when the dimrization 
is parametrized by the interaction for the alternating spin chain (1,|). 
The shift in the Schottky-like peak is to higher T, in contrast to the shift to 
lower T in Fig. (3). A large dimerization pushes the AF peak much farther 
than the F peak, resulting in a two-peak structure. 

Figure 6: The heat capacity for the alternating spin chain (§,§) with 
the parametrization as in Fig. (5). The F and AF peaks are now very 
dramatically separated. 

Figure 7: The xT vs T in ferrimagnetic chains: (a) for the system (1,^), 
(b) for (|, |) and (c) for (§, 1), for various values of the dimerization parame- 
ter 5. That dimerization enhances the AF character of ferrimagnets is clearly 
observable here. 

Figure 8: The free energy for the alternating spin chain (1,|) as a function 
of temperature for several values of the frustration parameter a. The same 
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behavior was found for (|, |) and (|, 1). 

Figure 9: The dependence of the excitation energies on frustration pa- 
rameter a for the spin system (1,|) shows the mode softening that heralds 
the transition to spiral state. 

Figure 10: Specific heat vs T : (a) for the spin system (1,|), (b) for (|, |) 
and c) for (|, 1), for various values of the frustration parameter a. 

Figure 11: The behavior of \T vs T of the ferrimagnetic system (1,|) on 
a chain for different values of a. 
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